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We study the heat capacity of isolated giant vortex states, 
which are good angular momentum (L) states, in a meso- 
scopic superconducting disk using the Ginzburg-Landau (GL) 
theory. At small magnetic fields the L=0 state qualitatively 
behaves like the bulk sample characterized by a discontinuity 
in heat capacity at T^. As the field is increased the disconti- 
nuity slowly turns into a continuous change which is a finite 
size effect. The higher L states show a continuous change 
in heat capacity at Tc at all fields. We also show that for 
these higher L states, the behavior of the peak position with 
change in field is related to the paramagnetic Meissner effect 
(irreversible) and can lead to an unambiguous observation of 
positive magnetization in mesoscopic superconductors. 

PACS numbers: 74.25.Ha; 74.60.Ec; 74.80.-g 



Superconducting state in mesoscopic samples can have 
strikingly different properties as compared to the corre- 
sponding bulk samples because by definition [1-10], the 
intrinsic length scales in them, like the penetration depth 
and the coherence length are of the order of the sample 
dimensions. Essentially, in a magnetic field a bulk type 
II superconductor can exist either in the Meissner state 
(macroscopic condensation of the Cooper pairs to the 
zero momentum state) or it can exist in a mixed state 
of higher momenta with infinite components (when there 
are vortices in the sample). In the presence of a bound- 
ary there can be another state called the giant vortex 
state wherein a very thin region near the surface is super- 
conducting ||ll|. This edge state in bulk samples affects 
resistivity measurements because contacts applied at the 
edges of the superconductor are shorted by it. There- 
fore, resistivity measurements are useful to determine the 
magnetic field to nucleate the giant vortex states. Oth- 
erwise these edge states are not capable of affecting the 
thermodynamic properties of the sample inside the phase 
boundary [|2j. Besides, in a large sample it is not possi- 
ble to isolate the individual giant vortex states because 
the sample encloses a large amount of flux and a large 
number of different giant vortex states are equally well 
formed at the same applied field. 

Recent developments in nano-fabrication techniques 
now allow us to reach the limit of isolating these indi- 
vidual giant vortex states in mesoscopic samples. We 
expect the thermodynamic properties of these isolated 
giant vortex states (that can be designated an angular 
momentum quantum number L) to be completely dif- 
ferent from each other and also from the bulk state. Re- 
cently, the phase boundary and magnetization |^-||] 
of these isolated L states have received some attention. 



Magnetization in superconductors is due to a gauge de- 
pendent equilibrium current and thus it is different in 
origin from magnetization of magnetic materials where 
the ordering of individual dipole moments, their dynam- 
ics and their internal energy scale determine the mag- 
netization. Heat capacity in superconductors is deter- 
mined by the thermal excitations of quasiparticles and 
therefore, it can eventually lead to a better understand- 
ing of the microscopic dynamics. For the first time we 
present a quantitative study of the heat capacity of the 
isolated giant vortex states, a quantity which is accessi- 
ble in an experiment. We show that at small fields there 
is a discontinuous jump in the heat capacity for L—0 
state at the transition temperature, a typical feature of 
second order phase transition, but at higher fields this 
turns into a continuous change. The line shape shows 
drastic changes during this crossover which can be inter- 
preted as a crossover from bulk-like to true mesoscopic 
behavior. On the other hand the higher L states show 
no discontinuity in heat capacity at all fields. Quanti- 
tative calculations strongly suggest that these features 
can be observed with present day experimental set-ups. 
Besides, we show that the paramagnetic Meissner effect 
(irreversible) pO| can be indirectly but unambiguously 
observed through the heat capacity measurements. 

In bulk samples Ginzburg-Landau (GL) theory is valid 
only close to the normal-superconductor transition when 
the gap at the Fermi surface approaches zero. So the 
jump in the specific heat of bulk samples at Tc agrees 
with the GL theory with high accuracy for some metals 
im . Well below the transition the observed specific heat 
scales, however, exponentially with temperature, signi- 
fying a unique energy scale associated with excitations 
and dynamics of quasiparticles in the system that orig- 
inates from the presence of a gap at the Fermi surface, 
and this does not come out in the GL theory. In the 
giant vortex state, the GL theory gives the correct mag- 
netization down to 0.33Tc and very small fields (7|j|]. The 
giant vortex states are current carrying edge states be- 
cause of which they are sometimes referred to as "gapless 
superconducting states" (for details we refer to ||ll|). In 
the absence of a gap the heat capacity will have a power 
law behavior below Tc according to the predictions made 
here based on the GL theory. The absence of gap can 
also be the reason why the GL theory correctly describes 
the magnetization at low temperatures p3| . 

We consider superconducting disks with radius R and 
thickness d immersed in an insulating medium. Taking 
the axis of the disk to be parallel to the z direction and 
keeping in mind the boundary condition at z = ±d/2. 
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one can expand the order parameter as a Fourier series 
'il){z,r) = T,kCos{^^)tpk{r)- Unlike in bulk samples is 
coordinate dependent in the x-y plane. Substituting in 
the GL equations one can check that for (Tr^/d)^ >> 1, 
the dominant contribution comes from '0fc=Oj which nat- 
urally means that if the thickness of the disk is less than 
the coherence length then the order parameter cannot 
vary in the z direction. Therefore we may assume a uni- 
form order parameter in the z direction and as a con- 
sequence the order parameter and currents become two 
dimensional with an effective penetration length that in- 
creases with thickness as Ae// = X^/d The flux ex- 
pulsion from the disk becomes negligible when Ae/ f ^ R- 
In that case the system is quantitatively described by the 
GL equation ||l5|, which is generally accepted as a 
Schrodinger like quantum mechanical equation 



(-iV - A)H = (1 - r)*(i - i^p), 



(1) 



that correctly describes the center of mass of the Cooper 
pairs. Here T is temperature in units of T'c(O) of bulk 
samples. Distance is measured in units of ^(0), the coher- 
ence length at zero temperature for bulk samples. The 
order parameter is measured in (T), the finite tem- 
perature order parameter of the bulk sample. The vector 
potential is measured in c?i/2e^(0) and the magnetic field 
in Hc2 = c?i/e£^(0). Our choice of units is the same as 
that in Ref . ||l6| . On the disk surface we require that the 
normal component of the current density is zero, which 
gives 



{-iV - A)n-^ - 0. 



(2) 



Rescaling the lengths (i? R\J{1 - T)) and fields {H 
H/ (1 — T)) we rewrite equation (1) in the following form: 
(-iV - AY"^ = ^{l - For the bulk samples the 

gradient term can be neglected jTll and in mesoscopic 
samples close to the normal superconductor transition 
the I^Pl^ term can be neglected [0J^]. We show that in 
the regime of our interest both the terms can be solved 
exactly in a semi-analytical way. 

The difference of the free energy G between the 
superconducting and the normal state, measured in 
Hl{Q)V/9>i:, can be expressed through the integral 



G 



2(1 -T){A- Ao).J- (1 - Tfm^] dr/V, (3) 



temperatures. As explained before for d < ^ equation (1) 
becomes two dimensional ]l4], i.e.. 



(4) 



We write the general solution as a truncated series in 
the solutions of the Linearized Ginzburg-Landau (LGL) 
equation (which is basically omitting the term in (4)) 



i.e., ^ = Ay^ro(f) -t- A'^"TL{r)exp{iL9). Here Voir) 
and rx,(r) are solutions of the LGL equation with eigen 
energies Ao and Xl, L being any non-zero integer giving 
the angular momentum quantum number. Hence, 



,1/2 



AL = (l + 2i.)-^-l, 



(5) 



where v is to be determined from the solution of the 
following equation. 

(L - ^)F{-u, ^ + 1' ^) - TTi^^'" + 1' ^ + = 

(6) 

where (p = HR^ is the fiux through the disk in units of 
fiux quantum hc/e. F{a,c,y) is the Kummer function, 
and finally 



TL{r)^r^expi- — )Fi^iy,L+l,—). (7) 

Although we retain only two terms in the series we will 
show that this is not an approximation in the regime of 
our interest. Now substituting ^ in equation (4) and 
simplifying we obtain 



XlAl 



a22A\ 



g^AqAl 



(8) 



(9) 



where Van = J T^dV, Va22 = J ^IdV and Vai2 = 
2 jTlrj^dV. We have to solve equations (8) and (9) to 
evaluate and Al- Three possible solutions are 



{Ao = Ao/aii, Al = 0), 



(^0 = 0, Al = XL/a22) 



(10) 



(11) 



over the disk volume V — nR^d, where Aq = Hr/2e^ 
is the external vector potential, and j — (5'*V5' — 
^V^*)/2j — I^PpA is the dimensionless supercurrent. 
The heat capacity difference between the superconduct- 
ing and normal states in units of Co = Hc{0)'^V/ (SttTc) is 
defined as C = -T^. The bulk limit of equation (3) is 



G = Hc{TyV/8TT Q which yields a specific heat jump 
of 2 in units of Co; whereas it can be accurately mea- 
sured to an order of magnitude smaller values at lower 



and 



(^0 



Aoa22 ~ ALai2 



011022 



a?2 



^ A^an^Aoai2 ^ (12) 

011022 — 0,12 



The solution in (10) corresponds to L=Q, the solution 
in (11) corresponds to L 7^ giant vortex states and 
the solution in (12) corresponds to a general mixed state 
solution. In (10) and (11), an contains all the correc- 
tion due to the non-linear term, Xi being eigenvalue of 



2 



the LGL equation. The mixed state solutions were dis- 
cussed earher |l7| . As we increase the number of terms 
in the expansion for '5 the solution (12) becomes more 
and more complicated and realistic but solutions (10) 
and (11) remain unchanged and exact. There is a defi- 
nite parameter regime where the giant vortex states have 
a lower free energy than the mixed vortex states and in 
this regime equations (10) and (11) give the same re- 
sults as the general numerical solutions |l^]. The stabil- 
ity analysis of these giant vortex states has been done 
analytically from the GL differential equations jl^] and 
numerically from the saddle point analysis of GL free 
energy [|l^, essentially proving their stability. Although 
the basic predictions of our work can be tested with just 
L=0 and the L—1 states which are under all conditions 
completely symmetric states the higher angular momen- 
tum giant vortex states have been detected experimen- 
tally jl^]. Also numerical calculations with a larger 
number of terms in the expansion, substituted into the 
GL free energy expression showed that numerically min- 
imized free energy corresponds to only one term in the 
expansion and the rest are identically zero in certain pa- 
rameter regimes. This parameter phase diagram is given 
in Ref. [ pO| and is the regime of our interest where the 
solutionsflO) and (11) are valid. 

In Fig. 1 (a), (b), and (c) we plot C/Cq vs T/TdO) 
at (p/(j)Q=0.1, 1.0 and 1.5 respectively, for all possible L 
states for a disk of i? = 4.0^(0). Figs. 1 (b) and (c) 
are displaced by 3 and 6, respectively, in the y-direction. 
The L—0 state is the ground state (shown in Fig. 2) but 
in decreasing fields it is possible to trap the system in 
the higher L states down to almost zero field due to the 
Bean-Livingston barrier pH , when fluxoids are trapped 
at the center of the sample and it can carry a large current 
1^,^. In Fig. 1 (a) the L—0 state shows a discontinuity 
at Tc (we refer to this as bulk- like behavior), while the 
higher L states show a continuous change at their corre- 
sponding TcS. As the flux is increased in Fig. 1 (b) and 
(c), the discontinuity in the heat capacity for the L—0 
state slowly changes to a continuous change, although 
there is a striking difference in the line shape of the L=0 
state from those of the L ^0 states. For high enough 
fields the line shape of the L=0 state will be the same as 
that of the L ^ states as will be shown in Fig. 3. A 
discontinuity in heat capacity is a characteristic feature 
of second order phase transition that can be observed 
in bulk samples. In finite samples the discontinuity is 
replaced by a continuous change. Symmetry breaking 
transitions, like solid-liquid melting transition, that are, 
however, first order transitions associated with a diver- 
gence in specific heat in the bulk, give a Gaussian curve 
in clusters [|3|. So the L—0 state in Fig. 1 (b) and (c) 
exhibits mesoscopic effect but not in Fig. 1 (a). In the 
L=0 state the Cooper pairs condense into a zero mo- 
mentum state {L being its angular momentum) just as 
bulk superconductivity is macroscopic condensation to 
the zero momentum state. Then in the absence of fields 
the boundary condition (2) for the currents at the bound- 



ary becomes irrelevant. Thus we get bulk-like behavior 
which persists at small but finite fields even though for 
infinitesimally small fields the effect of the boundary con- 
ditions set in. The higher L states show no discontinuity 
in the heat capacity at all fields. This line shape for the 
L=0 state in Figs. 1 (b) and (c) is intermediate between 
mesoscopic and macroscopic limits and would be an in- 
teresting feature to observe experimentally. 

In Fig. 2 we plot G/Gq (Go = 0AH^V/{8tt) vs (t>/(j)Q 
for a disk of i? = 4.0^(T) in dotted lines for different L 
values. The solid lines show {C/Co)peak, the magnitude 
of the peak value in C/Cq curves at corresponding fields 
for the different L states. The dashed curve would be 
{C /Co)peak for a bulk sample at zero field. The solid 
curves for any L appear approximately as a mirror re- 
flection of the dotted curves at the x-axis. The propor- 
tions of the solid and dashed curves are different in the 
y-direction because of which the crossings between the L 
states in the solid curves occur at the slightly different 
fields compared to that in the dotted curves. But the 
minimum in the dotted curve for a particular L is ex- 
actly at the same field of the maximum in the solid curve 
for the same L. This can be also argued from equation 
(3). Normally magnetic field in a bulk superconductor 
decreases the Cooper pair density and increases the free 
energy due to loss in the condensation energy associated 
with the breaking of Cooper pairs. But mesoscopic sam- 
ples have a regime for each L state where the free energy 
decreases with increase in field (this is nothing but para- 
magnetic Meissner effect ||l^ because magnetization is 
related to the flux derivative of free energy) as can be 
seen from the dotted curves in Fig. 2. In this regime the 
Cooper pair density (l^'P) also increase with increase in 
magnetic field unlike that in bulk samples. For exam- 
ple in the inset to Fig. 2, we consider the L=2 state 
in a regime where it shows the paramagnetic Meissner 
effect. The solid curve gives the Cooper pair density pro- 
file along a line starting from the center of the disk and 
ending at the boundary of it at an applied fiux of 2(f>Q. 
The dotted curve gives the same at 40o- Hence the den- 
sity at every point can be enhanced by increasing the 
field and as a consequence free energy decreases with in- 
crease in field. In this regime heat capacity as well as 
the peak value of the heat capacity, at a fixed temper- 
ature increases with increase in magnetic field. Each L 
state shows the paramagnetic Meissner effect at low fields 
when the field density in the disk exceeds the applied 
field. Field density decreases linearly with the applied 
field and at higher fields in the diamagnetic regime 
the field density inside the disk is less than the applied 
field. In this monotonous field dependence the magnetic 
energy is proportional to the square of the field expelled 
(being always positive) and hence the free energy mini- 
mizes when Cooper pair density maximizes and field den- 
sity is the same outside as well as inside the sample. In 
fact any observable quantity (like heat capacity) that de- 
pends on the density can identify the minima of the free 
energy and hence the regime where it decreases with in- 
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crease in field or the regime of paramagnetic Meissner 
effect. Hence at any fixed temperature, if the heat ca- 
pacity of a mesoscopic sample shows enhancement with 
increasing field, it can be taken as an indirect but unam- 
biguous observation of the paramagnetic Meissner effect 
(irreversible). Such an unambiguous observation is not 
possible by a direct magnetization measurement as has 
been clearly shown in Fig. 11 in Ref. while the sam- 
ple is always diamagnetic the observed magnetization is 
most of the time paramagnetic. The solid curve for the 
L=0 state at small flux shows a sudden jump that is re- 
lated to the crossover from bulk-like to true mesoscopic 
behavior. 

For the higher L states the difference in {C /Co)peak 
in Fig. 2 becomes small for consecutive L states but this 
will not create accuracy problems in measuring (C /Co) or 
{C/Co)peak of isolated L states up to considerably large 
L values, because the transition temperatures of the dif- 
ferent L states are different 0. Therefore, the (C/Cq) 
for consecutive L states peak at different temperatures 
and there can be a large difference between {C/Cq) as 
well as {C'/Co)peak of consecutive L states at particular 
temperatures which will help to isolate them experimen- 
tally. For example, we show the (C/Cq) for all possible L 
states at 4>/(f>o—llA in Fig. 3. The ground state at this 
field is the L—3 state (shown in Fig. 2). The difference in 
(C/Co) for consecutive L states at certain temperatures 
is accentuated by arrows and it can be much larger than 
the corresponding differences in (C/Co)peafe- The figure 
also shows that at large fields there is no noticeable line 
shape difference between the L—0 and L giant vortex 
states. 

Metastability can create some problems in an exper- 
imental situation but one faces the same problem in 
magnetization experiments. However, by systematically 
sweeping up and down the field one can perform measure- 
ments in almost the entire relevant part of phase space. 
One can just measure heat capacity as a function of mag- 
netic field and in turn obtain information about heat ca- 
pacity as a function of temperature. Ref. [Tc| ] also proves 
that by varying magnetic field at constant temperature 
or by varying temperature at a constant field one can 
arrive at the same stable state. So keeping the experi- 
mental set up unchanged one has to replace the detector 
by a calorimeter. 

In summary, we have shown that superconducting 
disks can give us an opportunity to observe a crossover 
of a second order phase transition in bulk system to such 
a transition in a finite system, where the discontinuity 
in heat capacity at the transition temperature is grad- 
ually smoothed out, accompanied by drastic changes in 
the line shape. The giant vortex states with non-zero 
angular momentum can exhibit finite size effects like a 
continuous change in the heat capacity at the transition 
temperature. Heat capacity measurements can also lead 
to an unambiguous verification of the existence of param- 
agnetic Meissner effect (irreversible) in mesoscopic sam- 
ples. Besides, we show that the GL theory can be suf- 



ficiently simplified to incorporate mesoscopic corrections 
to heat capacity because it retains accurately the spatial 
variations of order parameter that is determined by the 
coherence length and quiet independent of the penetra- 
tion length for Ag/ / > i?. 
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Figure captions 
Fig. 1 (a) C/Co versus T/Tc(0) for a disk of radius R = 
4^(0) at a flux 0/0o=O.l for all possible L states, (b) 
The same as (a) but (p/(f)Q=1.0. (c) (j)/(po=^-5. 
Fig. 2. The solid curves show {C/Co)peak versus (fy/do for 
all possible L states. The dotted curves show (G/Gq) ver- 
sus 4>/(l)o for all possible L states. The first few L states 
are labeled with their corresponding L values. Radius R 
of the disk=4.0^(r). The dashed curve gives {G/Go)peak 
for a bulk sample at zero field having the same coherence 
length. The inset shows the density profile along a radial 
line of the disk at two different fields. Co and Go, de- 
fined in text, contains all sample parameters. Thus the 
figure holds for all samples and at any temperature, the 
radius being four times the coherence length at that tem- 
perature. The radius in microns has to be thus properly 
chosen using the equations defined immediately below 
equation 2. 

Fig. 3. C/Go versus T/Tc(0) for a disk of radius R = 
4^(0) at a flux (p/cpo^llA for all possible L states. 
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